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I .  INTRODUCTION 


The  free  electron  laser  converts  the  kinetic  energy  of  a  relativistic 
electron  beam  into  a  coherent  electromagnetic  field  through  the  combined 
action  of  a  transverse  wiggler  field  and  the  electromagnetic  field.  For 
scaling  the  free  electron  laser  to  very  high  average  powers,  it  is  necessary 
to  operate  the  laser  in  a  Master  Oscillator  -  Power  Amplifier  (MOPA) 
configuration.  The  two  most  important  sources  of  high  energy  electron  beams 
suitable  for  use  in  the  FEL  are  the  rf  linear  accelerator  and  the  induction 
linear  accelerator.  The  rf  linear  accelerator  technology  is  mature  and  the 
linac  is  capable  of  producing  beams  with  peak  currents  of  a  few  hundred 
amperes  sufficient  for  sustaining  oscillation  in  a  resonator.  In  addition, 
the  electron  beam  brightness  from  the  rf  linear  accelerator  is  ct  least  an 
order  of  magnitude  greater  than  what  is  required  for  successful  FEL  operation. 
Moreover,  the  optical  output  obtained  from  an  rf  linac  FEL  consists  of 
extremely  short  micropulses  that  enable  the  laser  beam  to  propagate  through 
the  atmosphere  without  stimulated  rotational  Raman  scattering  (SRRS)  by 
This  particular  feature  has  important  implications  for  the  BMD  missions  for 
which  the  free  electron  lasers  are  being  seriously  considered.  However,  for 
successful  amplifier  operation  with  sufficiently  high  gain  (>15  dB)  and  high 
extraction  efficiency  (>20%),  it  is  desirable  that  the  peak  current  of  the 
accelerator  be  in  the  1  kAmp  or  greater  range. 

In  this  report,  we  discuss  a  potential  method  to  increase  the  peak  current 
of  a  rf  linear  accelerator  without  significantly  degrading  the  transverse 
emittance  of  the  beam.  This  can  be  achieved  at  the  expense  of  increasing  the 
longitudinal  phase  space  of  the  beam;  if  this  increase  in  longitudinal  phase 
space  is  acceptable,  the  FEL  performance  does  not  degrade  much.  In  the  next 
section  we  describe  our  concept  for  increasing  the  micropulse  peak  current. 
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We  then  estimate  the  change  in  brightness  of  the  beam  using  first  order 
transport  properties  of  the  beam.  In  Section  III,  we  anlayze  the  performance 
of  the  PEL  amplifier  for  different  electron  and  optical  beam  parameters. 
Section  IV  describes  our  results  of  transient  stimulated  Raman  scattering 
analysis  and  its  effect  on  both  up  and  down  links.  Section  V  summarizes  our 
findings  and  makes  recommendations  for  future  work. 
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II.  ANALYSIS  OF  BEAM  COMBINER 


Current  rf  linac  technology  can  provide  micropulse  peak  currents  of  up  to 

11  2 

500  amperes  at  a  brightness  of  >  10  amp/(m-rad)  .  For  amplifier  operation, 
a  threshold  current  of  -2  kamp  at  a  brightness  of  10  amp/(m-rad)  is 
typically  required.  At  higher  brightnesses  interesting  FEL  performances  can 
be  achieved  at  lower  threshold  currents.  In  what  follows,  we  shall  illustrate 
how  two  micropulses  can  be  combined  to  give  twice  the  current  at  twice  or 
slightly  less  than  twice  the  brightness.  The  concept  can  be  extended  in  a 
straightforward  manner  to  combining  more  than  two  pulses. 

There  are  several  schemes  for  combining  the  micropulses.  The  easiest  way 
would  be  to  combine  the  micropulses  from  several  near  identical  linacs  using 
"CHICANE".  This  technique  is  well  known  to  accelerator  physicists.  What  this 
requires  is  slight  energy  difference  between  the  different  accelerator 
micropulses.  The  CHICANE  works  essentially  as  an  electron  spectrometer  in 
reverse.  In  the  electron  spectrometer,  particles  of  different  kinetic 
energies  are  bent  through  different  angles  in  a  transverse  magnetic  field  and 
thus  can  be  separated.  If  the  particle  motions  are  reversed,  electrons 
starting  from  different  spatial  positions  with  differing  energies  can  be 
combined  to  form  a  single  pulse.  The  price  that  one  pays  for  increasing  the 
current  in  the  above  manner  is  that  the  energy  spread  of  the  beam  is 
increased.  If  the  increased  energy  spread  is  acceptable  to  the  FEL,  the 
concept  has  some  merit.  If  one  wants  to  reduce  the  accelerator  hardware  and 
use  only  a  single  accelerator,  one  could  combine  adjacent  micropulses  by  the 
following  technique.  At  the  end  of  the  linac  add  an  rf  cavity  operating  at  a 
frequency  such  that  the  successive  micropulses  enter  this  last  cavity  out  of 
phase.  Thus  if  the  first  micropulse  gains  energy  in  this  last  cavity,  the  2nd 
and  4th,  6th,  ...,  2Nth  pulses  will  lose  energy  while  all  odd  pulses  will  gain 
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energy.  This  is  illustrated  in  Fig.  2.1.  Clearly  if  the  frequency  of  this 

last  cavity  is  one  half  of  the  fundamental  frequency  of  the  linac,  one  can 

introduce  energy  spread  between  successive  micropulses. 

Other  frequencies  of  operation  for  this  linac  are  also  possible.  Let  f  be 

the  fundamental  cavity  frequency  of  the  accelerator  prior  to  the  final  cavity 

and  let  u  be  the  frequency  of  the  cavity  that  imparts  the  energy  dispersion  to 

the  different  micropulses.  If  the  harmonic  number  of  the  bunches  is  h,  then 

the  interpulse  duration  T  is  equal  to  h/f.  The  RF  phase  between  the 

p 

successive  micropulses  in  the  final  cavity  is  then  equal  to  2-kvI  .  Setting 

p 

this  equal  to  (2n+l)rr,  where  n  is  an  integer,  we  find  that 

w  -  (n  +  \  )f/h  (2.1) 

Choosing  n  to  be  equal  to  h  or  h-1,  we  get  a  frequency  for  the  final  cavity  to 
be  close  to  the  accelerator  frequency  f.  For  example,  if  f  -  500  MHz  and 
h  -  5 ,  we  see  that  u  -  500  ±50  MHz  would  make  every  other  micropulse  to  have 
a  slightly  larger  energy  than  the  ones  interleaving  them.  It  is,  of  course, 
not  necessary  to  choose  n  to  be  equal  to  h  or  h-1;  it  is  only  a  convenient  way 
of  sizing  the  cavity. 

Once  the  adjacent  micropulses  are  separated  in  energy,  it  is  a  simple 
matter  to  separate  them  spatially  with  a  dipole  magnetic  field.  The  two 
bunches  then  are  made  to  transverse  through  different  isochronous  and 
achromatic  transport  systems.  The  transport  systems  are  chosen  such  that  one 
bunch  travels  an  additional  distance  —  hjSc/f,  where  c  is  the  speed  of  light, 
and  f}  -  v/c=l.  These  two  transport  systems  bring  the  two  bunches  close  to 
each  other  to  positions  which  are  identical  to  the  ones  at  the  start  of  the 
transport  systems  (at  least  to  the  first  order).  The  two  bunches  then  can  be 
combined  to  form  a  single  bunch  by  reversing  the  process  which  separated  them 
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FIG.  2.1  -  SCHEMATIC  OF  THE  ACCELERATOR  CAVITIES  SHOWING  ENERGY  ADDITION 

TO  SELECTED  MICROPULSES.  BLACK  DOTS  REPRESENT  THE  POSITION  OF 
MICROPULSES  WITH  RESPECT  TO  THE  RF  WAVE. 


spatially  in  the  first  place.  The  schematic  of  one  such  beam  combiner  is 
shown  in  Fig.  2.2.  The  beam  combiner  therefore  can  be  separated  into  four 
components.  The  first  is  the  accelerator  cavity  tuned  to  the  frequency  v  to 
provide  energy  dispersion  between  adjacent  micropulses;  the  second  is  the 
magnetic  deflection  system  that  separates  adjacent  micropulses  spatially;  the 
third  is  the  isochronous  and  achromatic  transport  system  that  has  different 
path  lengths  to  bring  the  two  beams  to  coincidence  temporally;  the  last  is  the 
deflection  system  to  bring  the  two  bunches  together  spatially.  This  last 
component  is  the  mirror  image  of  the  second  component.  Once  the  adjacent 
beams  are  brought  together,  they  can  be  accelerated  further  in  cavities  tuned 
to  the  frequency  f.  This  will  reduce  the  effective  energy  spread  in  the  beam. 
For  example,  if  the  FEL  is  designed  to  operate  at  200  MeV,  then  the  beam 
combiner  can  be  made  to  operate  as  20±1  MeV.  The  2  MeV  energy  spread  imparted 
by  the  beam  combiner  amounts  to  only  12  of  the  final  energy  and  it  may  be 
within  the  energy  acceptance  of  the  tapered  wiggler  amplifier. 

For  the  magnetic  deflector  we  choose  a  system  first  proposed  and  analyzed 
by  Penner.1  He  analyzed  a  three-magnet  system  for  nondispersive  deflection  of 
a  nearly  parallel  beam.  Figure  2.3  shows  such  a  system.  The  curve  with 
arrows  represents  the  central  trajectory.  The  symmetry  axis  is  SS.  Electrons 
which  are  initially  parallel  at  the  entrance  phase  33  get  dispersed  at  SS  if 
they  are  of  different  energy  and  they  become  parallel  again  at  the  exit  plane 
00  so  if  we  break  the  deflection  system  at  SS,  we  should  be  able  to  spatially 
separate  the  two  adjacent  micropulses  which  have  different  energy.  Each  pulse 
then  is  transported  through  achromatic  and  isochronous  transports  to  the  S'  S' 
plane  where  the  beam  properties  are  identical  to  what  they  were  at  SS.  The 
difference  between  the  two  transport  systems  is  that  one  of  them  has  extra 
path  length  which  is  chosen  to  be  equal  to  the  inter-pulse  distance.  The 
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SCHEMATIC  OF  THE  TRANSPORT  OPTICS  FOR  THE  BEAM  COMBINER 
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THREE  MAGNET  SYSTEM  FOR  NONDISPERSIVE  DEFLECTION 


final  section  of  the  beam  combiner  is  the  same  as  the  component  between  SS  and 
00  in  Fig.  2.3. 

We  use  matrix  methods  for  analyzing  the  properties  of  the  magnet  systems. 

o 

Matrix  methods  for  calculating  trajectories  are  well  known  in  optics  and  have 
been  applied  to  electron  optics.  ’  '  In  the  matrix  theory,  the  position  and 
direction  of  travel  of  a  particle  entering  a  beam  line  is  represented  via  a 
vector  with  six  coordinates6 

x 
x' 

y 

y'  (2.2) 

i 
6 

The  coordinates  x  and  y  represent,  respectively,  the  horizontal  and 
vertical  displacements  at  the  position  of  the  particle,  and  x'  and  y' 
represent  the  angles  with  the  axis  of  the  beam  line  in  the  same  planes.  The 
quantity  2  represents  the  longitudinal  position  of  the  particle  relative  to  a 
particle  traveling  on  the  magnetic  axis  of  the  system  with  the  central 
momentum  designed  for  the  system.  The  remaining  quantity  5  -  (p  -  p  )/p 

O  0 

gives  the  fractional  deviation  of  the  momentum  of  the  particle  from  the 
central  design  momentum  of  the  system. 

The  effect  of  the  passage  of  a  particle  across  a  magnetic  element  or  a 
drift  space  is  represented  to  the  first  order  by  a  transfer  matrix  R.  The 
coordinates  X^  of  the  particle  at  the  end  of  the  element  are  then  given  in 
terms  of  those  at  the  beginning  X  as 

O 

X  -  R  X  (2.3) 

1  o 
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The  effect  of  successive  elements,  each  with  its  own  R  matrix  is  given  by  a 
total  R  matrix  R(t)  which  is  equal  to  the  product  of  the  individual  R 
matrices . 


R(t)  -  R  R  . . .  R  R 

n  n-1  2  1 


(2.4) 


The  transfer  matrix  R  is  a  6  x  6  matrix  and  the  problem  of  designing  the 
beam  combiner  reduces  to  (a)  finding  the  R  matrices  for  the  different  elements 
and  (b)  finding  a  suitable  coordination  of  the  R  matrices  that  give  the 
desired  total  R  matrix.  In  our  case,  we  would  like  R(t)  for  the  beam 
combining  to  be  as  close  to  a  unit  matrix  as  possible. 

Second  order  corrections  and  the  effect  of  deviation  from  the  nominal 
design  on  the  system  performance  may  then  be  applied.  If  all  the  bends  and 
focussing  are  carried  out  in  a  single  plane,  we  can  work  with  a  reduced  4x4 
R  matrix  to  make  a  preliminary  design  calculation.  This  will  be  the  approach 
that  we  take  in  this  report.  The  X  vector  in  this  case  will  have  only  four 
coordinates  x,  x' ,  i.  and  5. 

For  the  dispersive  section  53-SS,  the  transfer  matrix  R^  is  given  by 


0 

p/2sin/3 

0 

p(2cos/9-l) 

-2sin/3/p 

(4cos/3-3sec/8)/2 

0 

0 

2sin£(2cos/3-l) 

p(l-2cosa+2cos/?- 1 . 5sec/S) 

1 

p  (3,5-  2sina) 

0 

0 

0 

1 

where  a  is  the  magnet  bend  angle,  p  -  a/2  and  p  is  the  radius  of  curvature  of 
t.  e  central  trajectory.1 
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The  transfer  matrix  between  SS  and  S'S'  should  be  designed  to  be  unity. 
What  this  means  is  that  the  transport  system  be  achromatic  and  isochronous. 
Achromaticity  can  be  easily  achieved  by  demanding  that  the  total  phase  shift 
between  entrance  and  exit  be  an  integral  multiple  of  2x  radians.  A  little 
digression  may  be  necessary  to  explain  this  point.  A  particle  which  deviates 
from  the  central  trajectory  executes  betatron  oscillations  around  the  central 
trajectory.  To  first  order,  the  period  of  this  oscillation  depends  only  on 
the  central  momentum  and  magnet  parameters.  So  if  the  transport  system  is 
designed  such  that  the  total  length  equal  to  the  wavelength  of  the 
oscillation,  all  particles  will  end  up  with  the  same  x  and  x'  coordinates  as 
they  started  with.  Thus  the  transport  system  will  be  automaticaly  achromatic. 
However,  particles  with  different  energy  may  take  different  path  lengths 
through  the  transport.  So  the  system  may  not  be  isochronous.  Elements  with 
negative  dispersion  (momentum  compaction)  may  be  added  to  cancel  any  positive 
dispersion  present  in  the  system.  We  have  chosen  the  transport  between  SS 
and  S'S’  to  be  synthesized  from  elements  comprising  90°  deflection  system  as 
shown  in  Fig.  2.4.  The  transfer  matrices  for  this  system  can  be  written  as 


-1  2p(/2  -  1) 


n  p( 2  -  1/2) 


MAGNET  SYSTEM  FOR  NONDISPERSIVE  90  DEFLECTION 


'  -1  2pQ/2  -  1) 

0  *1 

-1/2  -p( 2  -  1/2) 

0  0 


0  0 

'0  0 

1  0 

0  1 


(2.7) 

where  the  superscripts  +  and  -  denote  whether  the  bend  is  to  the  right  or  left 
respectively  of  the  central  trajectory  at  the  entrance. 

The  transfer  matrix  for  path  A  can  be  written  as 


R+  R+  R*  R‘  R"  R+  R+  R'  R‘  R*  R+  R+  -  I 

dddddddddddd 


(2.8) 


where  I  is  the  identity  matrix.  Similarly,  the  transfer  matrix  for  path 
B  can  be  written  as 

R'  R+  R+  R'  -  I  (2.9) 

d  d  d  d 

The  validity  of  Eqs .  (2.8)  and  (2.9)  can  be  readily  seen  by  direct 

multiplication  using  Eqs.  (2.6)  and  (2.7).  If  Eq.  (2.9)  is  true,  it  is  easy 

to  show  that  Eq.  (2.8)  must  be  true  as  follows: 


R+  R+  R*  R'  R*  R+  R+  R'  R‘  R'  R+  R+ 
dddddddddddd 


R  R  R  R 

d  d  d  d 


R+  R*  R'  R+ 

d  d  d  d 


r'  r'  r+  r+ 

d  d  d  d 


R+  (R+  R'  R'  R+  )  R"  R'  (R+  R‘  R'  R+  )  R+ 

d  d  d  d  d  d  d  d  d  d  d  d 


-  R+  R'  R'  R+  -  I 

d  d  d  d 


(2.10) 
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The  trick  was  to  interchange  the  +  and  -  signs  in  the  middle  four  elements 
of  Eq.  (2.8)  since  they  multiply  out  to  identity  matrix.  The  rest  follows 

m 

naturally  by  regrouping.  If  we  denote  the  transfer  matrix  between  SS  and 
S' S'  as  R2,  we  have 

R,  “  I  (2.11) 

The  transfer  matrix  between  S'S'  and  00  can  be  easily  derived  to  be 

(2sina-l)2cosj8  p/ 2sin/3  0  p(2cos/J-l)  (l-2cosa)/2cos/J 

-2sin/?/p  0  0  2sirt/0(2cos/9-l) 

0  p(2cos/3-l)  1  p(3/3-2sina) 

0  0  0  1 


(2.12) 


The  total  matrix  R(t)  for  the  beam  combiner  system  therefore  is  given  by 


R(t)  -  R  R  R  -  R  R 

3  2  1  3  1 


(2.13) 


and  can  be  shown  to  be  given  by 


R(t)  - 


■1  p(2cot£-3coseca)  0 

0  1  0 

0  0  1 

0  0  0 


0 

0 

p(3a-4sina) 

1 


(2.14) 
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We  are  now  in  a  position  to  calculate  the  emittance  growth  in  the  beam 
combiner.  From  Eq.  (2.14),  we  have 

X  -  R  X  +  R  X'  and 

1  11  o  12  o 

X  '  -  X  '  (2.15) 

1  o  ' 

If  x  and  x  '  are  uncorrelated,  the  rms  emittance  is  given  by 

O  O 

2  2  1/2  2  2  1/2 
[  <x  >  <x'  >  1  before  the  beam  combiner  and  by  [  <x  >  <x'  >  1  after 
L  o  o  J  J  1  1  l  1 

the  beam  combiner.  We  have 

<x2>  <x/2>  -  <x2>  <x'2>  (R2  +  R2  <x'2>  /  <x2>  ]  (2.16) 

11  o  o  1  11  12  o  o 

The  emittance  growth  factor  therefore  is  given  by 

{R2  +  R  2  <x,2>  /  <x2>  }■  1/2 

-  {1  +  R2  <x'2>  /  <x2>  }  1/2  (2.17) 

Figure  2.5  shows  a  plot  of  the  emittance  growth  factor  as  a  function  of 
the  bend  angle  o.  In  plotting  Fig.  2.5,  we  have  assumed  that  the  radius  of 
curvature  to  be  25  cm  and  took  nominal  values  ofx  -1mm,  x'-l  mrad  and 

o  o 

Ap/p  -  12. 

The  normalized  electron  brightness  is  given  by 

21 

B  -  — - -  amp  /  mZ-  stered  (2.18) 

7  e  « 
x  y 

The  increase  in  pulse  length  is  given  by 

R34(62)1/2  (2.19) 
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25  cm,  ”0-1  mrad,  K0  -  t  mm,  (Ap/p)  -  0.01, 


EMXTTANCE  GROWTH  FACTOR  VERSUS  ALPHA 


The  peak  current  after  the  beam  combiner  is  given  by 


I  -  2  I  ct  /  (ct  +  R  <5Z>1/2) 

P  P  P  34 


(2.20) 


where  I  is  the  peak  current  in  each  micropulse  before  the  beam  combiner.  If 
we  assume  that  e  does  not  change  in  the  beam  combiner,  we  can  insert 
Eqs .  (2.20)  and  (2.17)  in  Eq.  (2.18)  to  obtain  the  brightness  improvement  as 


2c  t 

p 


1 


(ct 


R  <S2>1/2} 
34 


(R2  +  R2  <x'2>  /  <x2>  }1/2 

XI  12  o  'o 


(2.21) 


The  pulse  length  increase  is  plotted  in  Fig.  2.6  where  we  have  taken  the 
initial  pulse  length  to  be  -33  psec  (i  -  1cm).  Figure  2.7  shows  the 
brightness  improvement  as  a  function  of  the  bend  angle  a.  For  a  ~  60° ,  we  see 
an  improvement  of  80%  in  brightness. 

If  the  bend  angle  is  -  8°,  there  is  no  brightness  improvement.  For  bend 
angles  between  30  and  40° ,  the  brightness  improvement  is  -70%.  If  one  takes 
into  account  second  order  effects  and  any  misalignments  or  other  aberrations, 
this  increase  in  brightness  may  be  lost. 

We  have  shown  that  it  is  possible  to  increase  the  peak  current  and  increase 
or  maintain  the  beam  brightness  while  doing  this.  In  the  next  section  we 
discuss  how  this  affects  FEL  performance. 
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25  cm,  6 0  —  1  mrod,  XO  -  1  mm,  (ap/p)  "  0.01,  /“  1  cm 


LENGTH  INCREASE  VERSUS  ALPHA 


o 
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o 

LU 

a 


x 

a. 
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BRIGHTNESS  IMPROVEMENT  FACTOR  VERSUS  ALPHA 


III.  FEL  PERFORMANCE  ANALYSIS 


In  the  free  electron  laser,  the  electrons  are  given  a  periodic  transverse 
momentum  through  interaction  with  a  periodic  magnetic  bending  field  or  an  a 
electromagnetic  field  traveling  counter  to  the  motion  of  the  electrons.  Kroll 
and  McMullin7  have  shown  that  the  effect  of  these  two  fields  on  the  electron 
are  equivalent.  For  simplicity,  we  shall  consider  the  case  of  the 
relativistic  electron  beam  moving  inside  a  periodic  magnetic  field  structure 
called  the  wiggler  or  the  undulator.  The  magnetic  field  on  axis  can  be  either 
helical  or  planar  sinusoidal.  The  radiation  from  the  former  will  be 
circularly  polarized  while  in  the  latter,  it  will  be  linearly  polarized.  For 
simplicity,  we  shall  consider  planar  wiggler  whose  field  or  axis  is  specified 
by  the  vector  potential 


A 

W 


a  (z) 


dz' 


(3.1) 


The  radiation  field  propagating  in  the  z- direction  can  be  taken  to  be  derived 
from  the  vector  potential  A  given  by 


A 


a  (z)  sin  (k  z  -  w  t) 
s  s  s 


(3.2) 


The  quantities  a  and  a  are  the  dimensionless  vector  potentials.  For 

W  8 

magnetostatic  wigglers,  a  is  typically  of  the  order  of  unity,  while  for 

10  2 

radiation  at  a  wavelength  of  1  /^ra  and  a  flux  of  10  w/cm  ,  a  is  of  the  order 

-  4 

of  10  .  Thus,  for  practical  parameters  of  interest,  a  can  be  neglected 

compared  to  a  .  We  should  point  out  that  the  assumption  that  A  does  not 

depend  on  transverse  coordinates  cannot  be  strictly  satisfied,  but  if  the 
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transverse  dimensions  of  the  electron  beam  are  small,  Eq.  (3.1)  is  a 
reasonable  approximation  for  field  close  to  the  axis  of  the  wiggler.  Later, 
we  shall  consider  off-axis  effects  on  the  performance  of  the  free  electron 
laser . 

The  assumption  that  the  vector  potential  A  is  independent  of  transverse 
coordinates  implies  that  the  transverse  canonical  momentum  is  a  constant  of 
the  motion.  The  transverse  velocity  of  the  electron  can  therefore  be  written 
immediately  as 

a  .  a 

—  -  B  -  j/21  — —  sin  k  (z')  dz*  -  V?  — —  sin  (k  z  -  w  t)  (3.3) 

C  x  7  J  w  7  a  a 


where  we  have  assumed  that  the  electron  has  no  transverse  velocity  on  entrance 
into  the  wiggler. 

The  energy  equation  is  given  by 


<*7 

dt 


e 

me 


(0.E) 


(3.4) 


The  electric  field  E  can  be  obtained  from  Eq.  (3.2)  by  straightforward 
differentiation.  Using  Eq.  (3.3)  for  /S  we  can  write  Eq.  (3.4)  as 


d7 

dt 


w  a  a  r 

saw 

—  51 
f  w 


in  ij>  +  sin  -  2k  z  +  2w  t) 

S  S 


where 


a 

C 


(2  k  z  u  t) 

S  8 


+  k  z 

8 


w  t 

8 


(3.6) 
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The  angle  ip  denotes  the  phase  of  the  transverse  velocity  of  the  electron 


relative  to  the  radiation  field.  The  evolution  of  ip  is  obtained  by 

* 

differentiation  of  Eq.  (3.6)  as 


d0 

dz 


k  (z)  +  k 

**  * 


(3.7) 


If  sin  ip  in  Eq.  (3.5)  is  positive  and  quasi-stationary ,  the  first  term  on  the 
right  hand  side  of  Eq.  (3.5)  will  contribute  to  the  slowing  down  of  the 
electron.  The  second  and  third  terms  are  oscillating  at  about  half  the  period 
of  the  wiggler  and  do  not  contribute  to  any  significant  net  flow  of  energy 
between  the  electrons  and  the  radiation  field.  For  this  reason,  these  terms 
are  usually  neglected.  The  condition  when  f  -  0  is  called  the  resonance  or 
syncronism  condition.  If  we  substitute  w  -  k  c  in  Eq.  (3.7),  we  obtain 

k 

ip'  -  k  -  — —  (1  +■  a  2  -  2a  a  cos  ip)  (3.8) 

w  „  2  *»  w  s 

27 

2 

where  we  have  used  /3  -  1  whenever  permissable  and  (1  -  f)  )  -  1/27  •  Setting 

z  z  z 

V>'  -  0,  the  resonance  condition  is  obtained  as 


k  -  27Z  k  (z)/(l  +  a  Z) 

aw  w 


(3.9) 


Thus,  the  motion  of  the  electron  in  7  and  ip  space  is  determined  by  the 
solution  of  Eq.  (3.8)  and 


T' 


k  a  a 

s  s  »  . 

-  sm  ip. 

7 


(3.10) 


The  solution  of  Eqs .  (3.8)  and  (3.10)  for  different  initial  conditions  of 
and  7  ,  and  averaging  appropriately  .over  these  conditions  will  enable  us  to 

O 

determine  the  performance  of  the  free  electron  laser. 

Along  with  the  set  of  Eqs.  (3.8)  and  (3.10)  we  have  to  solve  the  equations 
for  the  evolution  of  the  optical  field.  If  we  let  k  -  k  +  5k  ,  it  is  easy 

s  so  s 

to  show  that 

2imeZ  \ 

5k  - - 

s  mew  a 

8  8 


N-4] 


(3.11) 


and 


a 


8 


2ffne2  fsin  i/> 

mew  w  I  7 
s  v 


(3.12) 


where  the  bar  denotes  averaging  over  all  the  electrons . 

Equations  (3.8),  (3.10),  (3.11)  and  (3.12)  are  solved  numerically  with 
k  (z)  and  a  (z)  being  known  input  functions  of  z.  For  N  electrons, 
Eqs.  (3.8)  and  (3.10)  really  constitute  a  total  of  2  N  equations. 

The  effect  of  emittance  on  the  performance  of  the  FEL  can  be  investigated 
by  means  of  an  effective  energy  spread.  The  energy  spread  is  taken  into 
account  by  taking  an  appropriate  distribution  of  the  initial  conditions  of  the 
electrons . 

The  FEL  has  many  physical  parameters  that  can  be  varied;  it  is,  therefore, 
impossible  to  make  a  complete  parametric  study  of  the  performance  of  the  FEL. 
To  focus  our  problem  on  the  effect  of  beam  combining  on  the  amplifier 
performance,  we  had  to  limit  the  independent  parameters  to  only  a  few.  A 
natural  limitation  is  the  length  of  the  wiggler.  Since  the  rf  linac 
micropulse  lengths  are  very  short  (<30  psec)  ,  long  wigglers  are  susceptible  to 
the  optical  slippage  problem.  The  electron  beam  slips  behind  the  optical  beam 
by  one  optical  wavelength  for  every  wiggler  period  traveled.  If  this  slippage 
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amounts  to  >10X  of  the  pulselength,  it  is  clear  that  the  electron  beam  tail 
would  not  contribute  much  to  the  generation  of  coherent  light.  Very  long 

wiggler  designs  are  based  on  the  principle  of  optical  guiding  (refractive 
guiding)  wherein  the  electron  beam  acts  like  a  fiber  and  guides  the  optical 
beam  through  the  wiggler  preventing  the  diffraction  expansion.  To  study  the 
physics  of  optical  guiding  correctly,  one  has  to  solve  the  paraxial  equation 
in  three -dimens ions .  Since  we  are  presenting  only  the  results  of  a 

one -dimensional  calculation,  we  did  not  want  to  invoke  the  optical  guiding  to 
increase  the  extraction  efficiency.  Accordingly,  we  restricted  the  wiggler 
length  to  three  times  the  Rayleigh  length  of  the  incident  optical  beam. 

The  second  restriction  that  we  have  placed  on  our  calculations  is  that  the 
wiggler  period  be  constant.  While  varying  the  taper  of  the  wiggler,  one  can 
change  A  ,  a  or  both.  Current  day  designs  have  all  opted  to  keep  A  constant 
because  of  the  ease  of  construction.  a^  is  varied  by  changing  the  gap  between 

the  pole  pieces.  In  what  follows,  we  have  taken  A  to  be  4  cm  and  the  optical 

wavelength  to  be  1  /um.  The  wiggler  length  is  taken  to  be  20  m. 

The  magnetic  vector  potential  (normalized  a  )  is  taken  to  vary  as 

a  -  a  exp[-(az)2]  (3.13) 

w  WO 

over  the  length  of  the  wiggler.  This  exponential  form  has  the  advantage  that 
the  wiggler  is  relatively  untapered  at  the  entrance  (i.e.,  z/L  «1)  while  the 
taper  is  appreciable  towards  the  end  of  the  wiggler.  This  allows  the  small 
signal  gain  to  build  up  very  fast  in  the  front  section  of  the  wiggler  as  well 

as  for  the  adiabatic  capture  of  as  many  electrons  as  possible  in  the 

ponderomotive  potential  well.  The  final  resonant  energy  of  the  electrons  is 
given  by 
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7R£  “  [Aw{X  +  \o  exp[-2(aL)2]}/2Aj1/Z  (3.14) 
while  the  resonant  energy  at  wiggler  entrance  is 

41  +  a«}/2A.]12  (3-15) 

The  energy  taper  therefore  can  be  written  as  (7  -  7  )/7  •  This  is  usually 
expressed  as  a  percentage.  The  smaller  the  value  of  a  in  Eq.  (3.13),  the 
smaller  is  the  energy  taper.  For  a  given  input  optical  flux,  the  intial 
trapping  fraction  increases  with  decreasing  a  while  the  deceleration 
efficiency  increases  with  increasing  a.  The  extraction  efficiency,  which  is 
approximately  the  product  of  trapping  fraction  and  the  deceleration  efficiency 
reaches  a  maximum  for  a  particular  value  of  a  that  depends  on  other  e-beam, 
wiggler  and  optical  parameters. 

The  two  most  important  parametric  studies  that  are  relevant  for  this 
concept  are  the  effects  of  peak  current  and  energy  spread  on  extraction 
efficiency  and  gain.  Figure  3.1  shows  a  plot  of  the  extraction  efficiency  as 
a  fraction  of  aL  for  three  different  peak  currents.  From  this  figure,  it  is 
clear  that  for  the  taper  function  that  we  have  chosen,  the  optimum  value  of  aL 
is  -0.8  for  peak  current  in  the  range  indicated.  Current  rf  linacs  are 
capable  of  providing  peak  current  of  400  to  500  amp  with  good  brightness.  For 
the  beam  combiner  concept  to  work  well,  we  would  like  the  peak  current  to  be 
in  the  kiloampere  range.  In  the  calculations  presented  above,  we  have  assumed 
that  the  incident  energy  distribution  is  gaussian  with  a  one  sigma  point  of 
0.5  percent.  This  represents  a  FWHM  energy  spread  of  approximately  one 
percent.  A  1  MeV  spread  imposed  by  the  beam  comnbiner  at  20  MeV  would 
correspond  to  one  percent  of  energy  spread  at  100  MeV,  which  is  close  to  the 
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EFFECT  OF  HIGGLER  TAPER  ON  EFFICIENCY 


nominal  energy  for  these  calculations.  Based  on  Fig.  3.1,  we  kept  aL  -  0.8 

while  studying  the  effects  of  other  parameters  on  the  performance  of  the  FEL 

*• 

amplifier.  Excursions  around  this  nominal  value  were  also  carried  out. 

Figure  3.2  shows  a  plot  of  the  efficiency  as  a  function  of  peak  current  for 
five  different  values  of  aL.  As  aL  is  increased,  the  current  necessary  to 
achieve  a  given  extraction  efficiency  increases  for  low  values  of 
efficiencies.  In  the  parameter  range  that  we  have  investigated,  we  find  that 
maximum  efficiency  is  obtained  for  aL  -  0.8  for  a  peak  current  of  1200 
amperes.  What  is  interesting  about  Fig.  3.2  is  that  the  shape  of  the  curves 
resembles  S ,  indicating  a  sudden  increase  in  the  extraction  efficiency  beyond 
a  threshold  value.  For  example,  take  the  case  of  aL  -  0.8.  Between  100  and 
500  amperes,  the  efficiency  increases  linearly  from  0.2  percent  to  1  percent. 
The  efficiency  jumps  from  a  value  of  -1  percent  at  500  amperes  to  -A  percent 
at  600  amperes  and  -7  percent  at  700  amp.  Beyond  700  amperes,  the  extraction 
efficiency  still  increases  but  the  increase  is  less  dramatic.  The  primary 
reason  for  attempting  to  increase  the  peak  current  is  evident  from  this  figure 
as  the  extraction  efficiency  increases  faster  than  linearly  beyond  the 
threshold. 

This  increase  in  extraction  efficiency,  of  course,  comes  at  the  expense  of 
increasing  the  longitudinal  phase  space.  It  is  therefore  of  some  interest  to 
find  what  the  penalties  are  of  increasing  the  energy  spread  in  the  beam. 
Figure  3.3  shows  a  plot  of  the  extraction  efficiency  as  a  function  of  energy 
spread  (FWHM) .  The  total  energy  spread  would  be  about  twice  of  what  is  shown 
in  the  figure.  The  efficiency  drops  fairly  linearly  from  a  maximum  for  zero 
energy  spread  to  almost  zero  for  a  value  of  (A7/7)  “  1.82.  As  the  peak 
current  is  increased,  the  value  of  (A7/7)  at  which  a  given  extraction 
efficiency  is  reached  increases  indicating  the  utility  as  increasing  the  peak 
current.  It  is  also  clear  from  the  figure  that  if  the  beam  combiner  operates 
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EFFECT  OF  ENERGY  SPREAD  ON  EFFICIENCY 


at  a  lower  energy  than  the  operating  energy  as  the  FEL,  the  net  percentage 
energy  spread  is  lower.  Conventional  linacs  operating  with  micropulses 
spanning  an  angle  6  of  the  rf  wave  have  an  energy  spread  “0  /8  due  to  the 
sinusoidal  variation  of  the  accelerating  field.  For  0  -  0.2  rad,  the  energy 
spread  would  be  -0.5  percent.  If  the  additional  energy  spread  introduced  by 
the  beam  combiner  is  of  the  same  order  as  this  or  less,  the  loss  in  efficiency 
due  to  increased  energy  spread  would  be  small  compared  to  the  increase  in 
efficiency  due  to  increased  peak  current. 

Finally,  the  effect  of  input  optical  flux  on  the  extraction  efficiency  was 

2 

studied.  Increasing  the  input  flux  from  1  to  10  Gw/cm  improves  the 

extraction  efficiency  dramatically  as  can  be  seen  from  Fig.  3.4.  Since  we 

have  taken  the  Rayleigh  length  to  be  one  third  of  the  wiggler  length, 

2 

a  peak  flux  of  10  Gw/cm  corresponds  to  a  power  of 
1010  x  ^  x  10  6  x  104  »  700  Mw.  If  the  duty  factor  of  the  linac  is  «10  3, 
this  would  require  an  oscillator  with  less  than  1  Mw  average  power.  The  peak 
output  power  would  be  in  the  range  of  10  -  30  Gw  giving  a  stage  gain  of  15  to 
100  (12  to  20  dB)  depending  on  the  input  power. 

Our  investigations  show  that  a  successful  beam  combiner  can  bring  immense 
benefits  to  the  operation  of  the  FEL  in  a  MOPA  configuration.  Longer  wigglers 
that  take  advantage  of  the  optical  guiding  can  lead  to  much  better  wiggler 
designs  and  improved  FEL  performance.  We  did  not  invoke  optical  guiding 
because  of  the  one -dimensional  nature  of  our  code.  Our  aim  was  to  show  that 
the  beam  combiner  has  many  advantages  in  the  operation  of  the  FEL.  Further 
analysis  remains  to  be  done  to  evaluate  the  MOPA  designs  and  their  performance 
thoroughly. 
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FIG.  3.4  -  EFFICIENCY  VERSUS  PEAK  CURRENT  FOR  DIFFERENT  VALUES  OF  INPUT  FLUX 


IV.  TRANSIENT  STIMULATED  RAMAN  SCATTERING 


In  the  downlink,  the  laser  beam  propagates  from  the  mission  mirror  to  the 
target  which  may  be  located  20  to  50  km  above  the  ground  at  the  time  it  is 
intercepted.  Atmospheric  absorption,  thermal  blooming  and  turbulence  are  not 
significant  issues  under  these  conditions.  However,  the  Stimulated  Rotational 
Raman  Scattering  (SRRS)  in  N2  in  the  atmosphere  imposes  important  constraints 
on  the  maximum  amount  of  flux  of  coherent  radiation  that  can  be  propagated 
through  the  atmosphere.  The  cross-section  for  the  SRRS  as  a  function  of  the 
wavelength  of  the  incident  radiation  has  been  measured  and  theoretically 
calculated  by  several  workers  .  Using  these  values  and  using  a  simple 
model  for  the  atmospheric  density  profiles,  it  is  possible  to  obtain 
steady- state  limits  on  the  laser  flux  that  can  be  transmitted  through  the 
atmosphere  both  in  the  uplink  and  in  the  downlink.  The  collisional  dephasing 
time  of  the  molecules  increases  with  increasing  altitude,  and  hence  for  pulsed 
lasers  the  Raman  scattering  during  downlink  is  usually  not  in  the  steady- state 
regime  but  rather  in  the  transient  regime.  The  transient  Raman  gain  is  less 
than  the  steady- state  Raman  gain  and  it  is  of  interest  to  find  how  the 
critical  altitude  (to  which  the  propagation  can  take  place  in  the  downlink 
without  significant  Stokes  conversion)  varies  as  a  function  of  the  pulse 
length,  wavelength,  laser  intensity  and  duty  factor.  In  this  section  we 
describe  how  we  have  solved  the  transient  Raman  gain  equations  numerically  for 
a  train  of  pulses  and  have  defined  the  critical  altitude  as  a  function  of  the 
average  laser  flux. 

The  pulse  formats  of  the  two  leading  free  electron  laser  concepts  are  shown 
in  Fig.  4.1.  The  individual  pulses  from  the  RFL/FEL  are  typically  20  to  30 
psec  long  and  are  very  short  compared  to  the  dephasing  time  in  the  molecules. 
However,  the  PRF  in  these  lasers  is  rather  high,  being  on  the  order  of 
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FIG.  4.1  -  TYPICAL  FEL  WAVEFORMS 
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100  MHz.  We,  therefore,  can  expect  the  molecules  to  have  memory  of  the 
passage  of  the  previous  pulses .  This  indeed  proves  to  be  the  case .  The 
quasi- steady  state  Raman  gain  for  a  train  of  pulses  is  higher  than  the  maximum 
Raman  gain  for  a  single  short  pulse  in  the  train.  The  1L/FEL  pulse  is 
typically  -50  nsec  long  and  is  somewhat  in  the  transient  regime.  However,  the 
interpulse  time  between  successive  IL/FEL  pulses  is  -100  psec,  which  is  long 
compared  to  the  collisional  dephasing  time.  Thus,  the  Raman  scattering  on 
each  pulse  takes  place  with  practically  no  memory  for  the  system  regarding  the 
propagation  of  earlier  pulses  through  the  medium. 

The  equations  governing  the  transient  stimulated  Raman  scattering  can  be 
written  as11 


( d/dz  +  a/viat)Ei(t,z)  -  i(2*wiVcZk1)NMfi*E.A  (4.1) 

( d/8z  +  3/v23t)E2(t,z)  -  i(2ww  Vc^  )NM  E  A*  (4.2) 

[d/at  +  rfi(z)]A*(t,z)  -  (i/ft)Mfi*(Pi  -  Pf)Ei*E2  (4.3) 

with  |  Ex  | 2  and  lEj2  being  the  intensities  of  the  pump  and  Stokes  fields, 

respectively;  A  is  the  material  excitation;  N  -  N(z)  is  the  number  density  of 

N2;  p ^  and  Pf  are  the  probability  densities  of  the  initial  and  final  states; 

M  -  <£  j  M |  i>  is  the  Raman  matrix  element  for  the  transition  between  |  i>  and 

|  f >  states;  and  r^(z)  is  the  dephasing  rate.  v^  and  are  the  group 

velocities  at  w  and  w  ,  and  we  can  assume  that  v  -  v  —  v. 

12  12 
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For  the  problem  that  we  are  trying  to  solve,  we  can  assume  that  the 


depletion  of  the  pump  radiation  as  well  as  the  induced  population  change  are 
negligible.  Then  the  above  equations  reduce  to 

(d/dz  +  3/v3 t)E  -  in  (z)E  (t  -  z/v)A*  ,  (4.4) 

2  11 

[3/at  +  r(z)]A*  -  -if?  E  *(t  -  z/v)E  ,  (4.5) 

2  1  2 


where 

f}i(z)  -  i(2frw22/c2^2)N(z)Mfi 
f7  -  M  ★(/>.  -  p  )/r 

2  f  i 

We  can  simplify  these  further  by  the  appropriate  change  of  variables:  let 
T  -  t-z/v  (retarded  time)  and  Z  -  z,  so  that  3/3Z  -  3/3z+3/v3T;  also,  let 
Ed . Z)  -  E2(T,Z),  F(T,Z)  -  ir,iEi(T)A*(T,Z),  H(z)  -  rj^Z)^,  and 

p(T)  -  (Ej2  0(T)  ,  where  Ei  is  a  constant,  and  fl(T)  is  the  step  function.  (It 
is  not  necessary  to  use  a  step  function  for  E^  in  our  procedure,  but  the 
computation  time  is  significantly  reduced  if  this  is  done,  and  the  equations 
are  somewhat  simpler;  for  the  problem  that  we  are  addressing,  the  actual  pulse 
shape  is  insignificant.)  Then  our  equations  can  be  simplified  to 

3E/3Z  -  F(T,Z)  (4.6) 

and 

3F/3T  -  P(T)H(Z)E(T,Z)-r(Z)F(T,Z)  (4.7) 
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There  exists  an  analytical  solution  of  these  equations  for  a  uniform 
11  12 

medium  '  ,  in  which  H  and  T  are  constant.  However,  in  the  atmosphere,  H  and 

T  are  functions  of  altitude  and  could  be  complicated  (especially  if  the 
temperature  variations  were  taken  into  account) ,  so  that  in  general  the 
analytical  approach  would  fail.  Several  numerical  procedures  were  considered 
and  rejected  by  us  as  inefficient.  Most  prominent  among  them  was  using  a 
conventional  finite  difference  scheme  on  a  two-dimensional  grid  to  solve  the 
equivalent  hyperbolic  equation.  This  method  would  require  producing,  storing, 
and  inverting  large  ill-conditioned  matrices,  and  if  even  at  one  point  there 
would  occur  a  numerical  overflow,  the  whole  inverted  matrix  would  be  incorrect 
(if  the  original  one  remained  nonsingular).  It  would  also  require  a  much 
finer  grid  (more  grid  points)  than  the  method  we  used.  Alternatively,  we 
could  use  a  numerical  two-dimensional  integration  procedure  over  rectangles  or 
hexagons  in  the  T-Z  plane  on  the  corresponding  double- integral  equation.  This 
method  would  require  too  rigid  and  inconvenient  relationships  between  the  step 
sizes  along  the  T  and  Z  directions. 

The  approach  that  we  took  was  based  on  the  Adam-Bashforth  type 

predictor -corrector  method,  which  is  conventionally  used  for  ordinary 

differential  equations.  We  have  chosen  the  fifth  order  predictor  and  sixth 

order  corrector  formulas  (we  could  ''afford"  such  a  high  order,  since  the 

Adams -Bashforth  type  methods  are  relatively  stable)  that  for  an  ODE 

13 

y'  -  f[x,y(x))]  have  the  following  form  : 

yn+i  -  yn+(h/720)(1901y'n-2774y'n  2-1274y'n_3+2511274y'n  /,)  (4.8) 

and 

y  -  y  +(h/1440)(475y'  +1427y'  -798y'  +482y'  -173y'  +27y'  )  (4.9) 

"n+l  -'ll  '  J  n+1  J  n  J  n-1  J  n-2  J  n-3  n-* 


36 


where  h  is  the  integration  step  size,  and  the  subscripts  refer  to  the  points 
on  the  x-axis  (separated  by  h) ,  at  which  y'  is  evaluated  and  y  is  estimated. 
The  main  advantage  of  using  a  high  order  predictor-corrector  method  is  in  its 
high  accuracy  at  low  cost  (since  h  can  be  chosen  to  be  relatively  large) .  The 
main  disadvantages  lie  in  the  numerical  instability  for  some  equations  and  the 
necessity  for  some  kind  of  single  step  method  to  initialize  the  multi: tep  one. 
The  accuracy  of  the  single  step  method  should  be  comparable  to  that  of  the 
multistep  one;  for  example,  the  Runge-Kutta  method  is  normally  used  to 
initialize  the  fourth  order  Adams-Moulton  multistep  one,  since  the  local  error 
in  both  of  them  is  proportional  to  h5.  (In  our  method  we  cannot  use  a 
Runge-Kutta  type  starter,  so  we  use  a  relatively  complicated  procedure, 
involving  the  Huen  method  and  several  increases  of  the  step  size.) 

We  start  our  calculations  from  the  following  initial  and  boundary 
conditions:  F(0,Z)  —  0,  E(0,Z)  —  E(T,0)  —  C,  where  C  is  some  positive 
constant  that  is  chosen  in  such  a  way  that  makes  a  numerical  overflow  or 
underflow  in  the  computations  unlikely.  Physically  this  means  that  we  are 
starting  the  SRRS  from  a  coherent  seed,  whose  actual  value  is  irrelevant  to 
us,  since  we  are  only  interested  in  the  Raman  gain,  which  equals 
in[E(T,Z)/E(T,0) ]  -  in[E(T, Z)/C] .  At  first,  F(T,0)  is  evaluated  at  the 
appropriate  points  from  the  analytical  solution  (which  is  easily  obtained, 
since  E(T,0)  is  a  constant).  Then  we  find  F  and  E  on  the  lines  (T,nAZ)  (n  is 
an  integer),  line  by  line,  until  we  span  the  region  of  interest  in  the  T-Z 
plane.  For  any  point  lying  well  inside  that  region  we  proceed  as  follows: 
using  the  predictor  formula  with  the  values  of  3F/3T  at  the  previous  five 
points  for  the  same  Z  (we  know  them  from  before),  we  "predict"  the  value  of  F; 
then,  using  the  corrector  formula,  with  the  value  of  F  that  we  just  found,  as 
well  as  the  values  at  the  previous  five  points  for  the  same  T,  we  find  the 
value  of  E  at  the  current  position.  Using  this  value  of  E,  we  evaluate  3F/3T, 


37 


and  use  it  in  the  corrector  formula  for  F;  we  repeat  this  process  twice. 
Thus,  we  predict  F  once,  correct  it  twice,  and  "correct"  E  three  times  (that 
is ,  use  the  corrector  formula) . 

When  the  pump  pulse  ends,  the  E  field  departs  from  the  region  of  interest 
with  the  speed  of  light,  while  the  material  excitation  (and  therefore  F) 
decays  due  to  collisions  between  the  air  molecules,  so  that  the  solution  to 
our  equations  becomes:  E(T+r,Z)  -  F(T,Z)exp(-T  r)  ,  where  7  is  the  duration  of 

C 

the  first  pulse,  r  -  T-T,  and  T  is  the  collisional  frequency.  Therefore,  the 
initial  conditions  on  F  for  the  next  pulse  will  depend  on  F(T)  and  time 
between  the  pulses.  This  variability  is  taken  into  account  in  our  program. 
If  a  very  long  pulse  needs  to  be  considered,  we  can  simply  set  the  time 
between  the  pulses  to  zero  (because  the  program  has  some  memory  limitations 
that  set  a  limit  on  the  single  pulse  duration) . 

The  results  of  the  program  are  the  values  of  the  total  and  the  local  gain 
G  (T,Z)  -  in[ E(T, Z)/'E(T,0) ]  and  G  (T,Z)  -  3(inE)/3E  -  F(T,Z)/E(T,Z) ,  which  can 

T  L 
later  be  plotted. 

In  order  to  test  the  performance  of  our  program,  we  evaluated  the 
analytical  solution  for  the  uniform  medium  case  and  compared  it  with  the 
corresponding  results  obtained  from  the  program.  The  agreement  was  to  at 
least  four  significant,  figures,  which  exceeds  the  resolution  of  most 
computer -generated  graphs.  We  also  set  the  interpulse  duration  to  zero  and 
ran  the  program  long  enough  to  reach  the  steady- state  limit  in  an  exponential 
atmosphere.  These  results  were  compared  with  analytical  results  that  one  can 
obtain  from  Eqs .  (4.6)  and  (4.7)  by  setting  the  time  derivative  in  Eq.  (4.7) 
to  zero.  Once  again  we  found  excellent  agreement  between  the  numerical 
results  and  analytical  solutions  verifying  the  stability  and  accuracy  of  the 
numerical  procedure. 
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The  transient  Raman  gain  equations  are  solved  to  yield  the  exponential  gain 

for  a  'square  hat'  pump  pulse  as  a  function  of  the  altitude  and  pump 

► 

intensity.  The  propagation  path  in  all  downlink  cases  is  assumed  to  be  from 
space  to  the  ground  at  an  angle  of  60“  to  the  zenith.  We  assume  that  an 

40 

intensity  gain  of  e  is  required  for  almost  complete  conversion  of  the  pump 

beam  into  the  Stokes  beam.  The  divergence  angle  of  the  Stokes  beam  is 

approximately  given  by  D/L  where  D  is  the  diameter  of  the  beam  and  L  is  an 

effective  gain  length  for  Raman  conversion.  Typical  values  for  D  and  L  might 

be  0.3  m  and  50  km,  respectively.  Thus,  the  converted  Stokes  radiation  is 

still  fairly  well  collimated.  However,  as  this  beam  propagates  downwards,  the 

gain  increases  because  of  increasing  N2  number  density  down  to  -45  km,  below 

which  point  the  gain  is  sensibly  constant.  The  conversion  from  the  first 

Stokes  radiation  into  the  second  and  higher  order  Stokes  beams  proceeds  very 

rapidly.  This  process  is  illustrated  in  Fig.  4.2.  Typical  conversion  lengths 

for  these  are  5  to  6  km.  The  higher  order  Stokes  beams  have  therefore  a  much 

larger  beam  divergence  angle  than  the  first  order  beam;  one  can  estimate  the 

distance  these  beams  would  travel  before  their  intensity  drops  by  a  factor  of 

-10.  To  take  this  physics  into  account  in  the  downlink  problem,  we  have 

40 

calculated  the  altitude  at  which  the  gain  becomes  e  and  then  subtracted 
10  km  from  this  altitude  to  obtain  the  critical  altitude  for  propagation  of 
the  laser  beam  with  a  given  intensity. 

The  calculations  were  carried  out  at  the  pump  wavelength  of  1.06  /an.  The 
stimulated  Raman  gain  scales  almost  linearly  with  the  pump  frequency  in  the 
wavelength  region  of  0.35  n m  to  2.0  /im.  It  is  therefore  a  simple  matter  to 
scale  the  intensity  linearly  with  the  wavelength  to  derive  the  results  for 
wavelengths  different  than  1.06  ^tm.  Figure  4.3  shows  the  critical  altitude 
down  to  which  the  RFL/FEL  pulses  can  propagate  as  a  function  of  the  duty 
factor  for  different  average  flux  levels.  The  micropulse  length  for  these 
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LIMITING  ALTITUDE  VERSUS  DUTY  FACTOR  FOR  RF  LINAC  FEL 


calculations  was  taken  to  be  30  psec  which  is  typical  for  an  accelerator 

working  around  500  MHz,  For  the  FELs  of  interest  here,  the  duty  factor  is 

expected  to  be  in  the  range  of  1  to  3  x  10  3 .  From  the  figure,  it  is  clear 

that  the  RFL/FEL  pulses  will  have  no  trouble  getting  down  to  -10  km  altitude 

with  significant  average  fluxes.  Figure  4.4  is  similar  to  Fig.  4.3,  except 

the  micropulse  length  is  set  at  100  psec.  In  the  latter  case  (i.e., 

Fig.  4.4),  the  Raman  gain  is  higher,  and  for  the  same  average  flux  and  duty 

factor,  the  critical  altitude  is  higher  than  that  for  the  case  of  shorter 

micropulse  length.  The  calculations  show  that  the  RFL  pulses  of  pulse  length 

less  the  100  psec  and  duty  factor  >  0 .  IX  will  have  no  trouble  propagating  down 

2 

to  -10  km  altitude  at  an  average  intensity  of  0.3  MW/cm  .  Figure  4.5  shows 
the  critical  altitude  as  a  function  of  peak  intensity  for  a  50  nsec  pulse  and 

CW  laser  beam.  A  typical  duty  factor  for  the  IL/FEL  is  -  5  x  10  .  Under 

these  conditions  the  minimum  altitude  down  to  which  the  pulse  will  propagate 
is  about  50  km  at  flux  levels  of  interest  to  SDI  missions.  There  are  two  ways 
of  bringing  this  altitude  lower.  The  first  one  is  to  operate  the  IL/FEL  at  a 
lower  peak  power.  The  penalty  one  must  pay  for  a  lower  peak  power  is  the 
lower  extraction  and  overall  efficiency  of  the  laser.  The  second  workaround 
for  the  downlink  SRRS  problem  is  to  use  optical  techniques  to  stretch  the 
pulse  length  of  the  output  beam,  thereby  reducing  the  peak  power.  Since  the 
50  nsec  pulse  has  a  lower  Raman  gain  than  a  true  CW  beam,  increasing  the  pulse 
length  makes  the  Raman  gain  approach  the  CW  gain  limit  more  and  more.  A  factor 
of  3  increase  in  pulse  length  does  not  change  the  results  appreciably.  To  get 
significant  decreases  in  the  critical  altitude,  the  pulse  length  has  to  be 
increased  by  at  least  a  factor  of  30. 

In  the  uplink,  the  laser  beam  goes  through  dense  air  in  the  begining  where 
the  relaxation  time  is  -  0.14  nsec.  In  this  case  the  steady- state  SRRS 
threshold  is  appropriate  for  the  IL/FEL  pulse  lengths  of  50  nsec,  while  the 
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UPLINK  SRRS  IRRADIANCE  LIMIT  VERSUS  DUTY  FACTOR 
(PULSE  LENGTH  EQUALS  30  PS) 


much  shorter  pulse  lengths  (e.g.,  =  30  psec)  of  the  RFL/FEL  can  experience 
transient  SRRS  gain  effects  on  the  uplink  propagation  path.  Thus,  transient 
SRRS  effects  must  be  considered  in  ' specifying  the  minimum  RFL/FEL  aperture 
size . 

Figure  4.6  shows  the  results  of  a  transient  SRRS  analysis  we  have  performed 
for  uplink  propagation  at  the  1.06  fim  wavelength  with  a  pulse  length  of 
30  psec,  which  is  typical  of  the  RFL/FEL.  The  peak  beam  irradiances  required 
to  achieve  uplink  SRRS  gains  of  20  and  40  nepers,  respectively,  are  plotted  as 
functions  of  the  duty  factor  in  Fig.  4.6  for  the  case  of  a  3  km  GBL  altitude 
and  40°  zenith  angle.  Operation  at  duty  factors  less  than  one  with  30  psec 
pulses  is  shown  to  rapidly  increase  the  threshold  for  significant  SRRS 
conversion.  (A  gain  of  40  nepers  corresponds  to  that  required  for  full 
conversion  of  the  uplink  beam  to  the  first- stokes  frequency,  with  the 
effective  beam  quality  being  seriously  degraded  as  the  beam  exits  from  the  top 
of  the  atmosphere  layer;  e.g.,  at  -  50  km  altitude.)  For  example,  considering 
a  20  neper  SRRS  gain  as  an  acceptable  limit,  allowing  for  the  presence  of 
sideband  energy  which  may  be  present  at  the  first  stokes  frequency  in  the 
RFL/FEL  output  spectrum,  the  threshold  irradiance  for  30  psec  pulses  is 
increased  to  about  20X  the  steady- state  threshold  value  for  duty  factors  of 
the  laser  less  than  10  2 .  Since  3  x  10  3  is  typical  for  the  RFL/FEL  duty 
factor,  based  on  present  technology  projections,  the  RFL/FEL  would  appear  to 
have  essentially  no  problems  with  uplink  SRRS  for  any  GBL  power  levels 
currently  of  interest. 
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WAVELENGTHS,  SCALE  PEAK  FLUX  WITH  X. 


LIMITING  ALTITUDE  VERSUS  PEAK  FLUX  FOR  LONG  PULSES 


V.  SUMMARY 


We  have  shown  in  this  report  that  the  FEL  performance  can  be  markedly 
improved  if  the  peak  current  can  be  increased  to  kiloampere  level.  These 
improvements  can  be  obtained  in  relatively  short  wigglers  (-20  to  25  m) 
without  utilizing  the  optical  guiding  for  increasing  the  extraction 
efficiency.  Therefore,  our  results  are  somewhat  conservative  in  this  sense. 
Gains  of  12  to  20  dB  can  also  be  realized  concurrent  with  extraction 
efficiencies  of  >  8X  in  the  amplifier. 

We  have  analyzed  on  scheme  for  combining  the  electron  micropulses  to 
increase  the  peak  current  without  degrading  the  beam  brightness.  The  analysis 
was  carried  out  to  first  order  using  matrix  methods.  The  preliminary  analysis 
was  limited  to  one  transverse  direction  only.  The  first  order  analysis 
clearly  indicates  the  feasibility  of  the  beam  combiner  concept.  It  requires 
further  verification  using  TRANSPORT  and  TURTLE  codes  that  are  available. 
Even  if  the  present  concept  is  found  to  have  drawbacks,  it  is  always  possible 
to  have  two  accelerators  that  operate  at,  say,  20  and  21  MeV  respectively  and 
combine  the  micropulses  from  these  two  accelerators.  The  resulting 
micropulses  can  be  further  accelerated  to  100  to  200  MeV  required  for  the  FEL 
MOPA  operating  in  the  visible  or  near  IR.  The  effect  of  aberrations  such  as 
misalignment  and  manufacturing  tolerance  on  the  performance  of  the  beam 
combiner  have  not  been  studied;  nor  have  we  included  second  and  higher  order 
effects.  Hopefully,  these  will  be  studied  in  a  future  investigation. 

We  have  also  investigated  the  down  link  transient  SRS  problem  through  the 
atmosphere.  The  RF  Linac  FEL  pulse  format  seems  to  be  ideal  in  reducing  the 
SRS  growth  rate  through  the  atmosphere.  For  flux  levels  of  interest  to  BMD, 
we  do  not  foresee  any  problem  in  propagating  the  optical  beam  through  the 
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atmosphere  on  the  targets.  In  the  uplink,  thermal  blooming  is  perceived  to  be 
the  driver  in  determining  the  size  of  the  optics  and  not  SRS  for  power  levels 
of  interest  in  the  BMD  missions. 
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